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We present a simple model of Texas hold'em poker tournaments which retains the two main 
aspects of the game: i. the minimal bet grows exponentially with time; ii. players have a finite 
probability to bet all their money. The distribution of the fortunes of players not yet eliminated is 
found to be independent of time during most of the tournament, and reproduces accurately data 
obtained from Internet tournaments and world championship events. This model also makes the 
connection between poker and the persistence problem widely studied in physics, as well as some 
recent physical models of biological evolution, and extreme value statistics. 



I. INTRODUCTION 

Physicists are now more then ever involved in the study 
of complex systems which do not belong to the tradi- 
tional realm of their science. Finance (options theory,...) 

human networks (Internet, airports,...) ||, the dy- 
namics of biological evolution ||, and in general of 
competitive "agents" || |[ 0] are just a few examples 
of problems recently addressed by statistical physicists. 
However, many of these systems are not isolated and are 
thus sometimes very difficult to describe quantitatively: 
a financial model cannot predict the occurrence of wars 
or natural disasters which certainly affect financial mar- 
kets, nor can it include the effect of all important exter- 
nal parameters (China's GDP growth, German exports, 
Google's profit...). Rather, these studies try to capture 
important qualitative features which, interestingly, are 
sometimes universal. In this context, universality means 
that large scale aspects of the real system are properly re- 
produced by a simple model which only retains the main 
relevant ingredients of the original physics. Adding fur- 
ther details to the model does not affect these universal 
properties. 

In the present work, we study a very human and play- 
ful activity: poker tournaments. Although a priori gov- 
erned by human laws (bluff, prudence, aggressiveness...), 
we shall find that some of their interesting properties 
can be quantitatively described. One of the appealing 
aspects of a poker tournament lies in the obvious fact 
that it is a truly isolated system, which is not affected 
by any external phenomenon. Two famous mathemati- 
cians (Emile Borel M, and later John von Neumann jj^j) 
contributed to the science of poker. However, they con- 
centrated on head-to-head games, like their most recent 
followers jig ], obtaining the best strategy in terms of 
the value of the hand and the pot. To our knowledge, 
the present work represents the first study of large scale 
poker tournaments. Note however that in a recent work 
||, the authors study head-to-head elimination tourna- 
ments involving seeded competitors, and apply success- 
fully their theory to the US college basketball national 
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championship. 

In the following, we introduce a simple model which 
can be treated analytically and which faithfully repro- 
duces some properties of Internet and live poker tourna- 
ments. Our main quantities of interest are the distribu- 
tion of the fortunes of surviving players, their decay rate, 
the number of different players owning the biggest for- 
tune at any given time during the tournament (dubbed 
the "chip leader"), and the distribution of their fortune. 
Interestingly, the constraint that a surviving player must 
keep a positive fortune relates poker tournaments to the 
problem of persistence O, [l^, |l3| , and the competitive 
nature of the game connects some of our results with 
recent models of competing agents || |], ||, |[ [?J . In ad- 
dition, the properties of the chip leader display extreme 
value statistics, a phenomenon observed in many physical 
systems § 0, @. 

In Section II, we define a stochastic model which re- 
tains the main identified ingredients of poker tourna- 
ments: i. the minimal bet grows exponentially with time; 
ii. players have a finite probability q to bet all their 
money. In Section III, we first solve the correspond- 
ing model for q = 0, which will allow us to make the 
connection with the persistence problem widely studied 
by physicists. In Section IV, we will show that q must 
physically take a specific value, and thus is not a free pa- 
rameter. The results of the model will compare favorably 
with actual data recorded from real Internet poker tour- 
naments and World Poker Tour main events. Finally, in 
the last Section V, we will consider the statistical prop- 
erties of the chip leader. In particular, we will show the 
connection with the "leader problem" arising in evolu- 
tionary biophysics, and the field of extreme value statis- 
tics which has recently attracted a lot of attention from 
physicists. 



II. A SIMPLE POKER MODEL 

Before addressing the basic rules of poker and the re- 
sulting definition of our model, we wish to introduce some 
useful poker terminology. In a real poker tournament, 
players first pay the same entry fee or "buy-in" (from 1 $ 
to 25000 $) which is converted in "chips" . Hence players 
are not betting actual money but chips. The total num- 
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ber of chips of a player is called his "stack" . At any time 
in the tournament, if a player decides to bet his entire 
stack, it is said that he is going "all-in" . 

We now describe the main aspects of a Texas hold'em 
poker tournament, currently the most popular form of 
poker. Initially, No players sit around tables accepting 
up to 9 = 10 players. In real poker tournaments, No 
typically lies in the range Nq ~ 10 — 10000. We do not 
detail the precise rules of Texas hold'em poker, as we shall 
see that their actual form is totally irrelevant provided 
that two crucial ingredients of the game are kept: 

• A tournament consists in a series of independent 
games or "deals". Before a deal starts, the two play- 
ers next to the dealer (i.e. the player dealing the cards) 
post the minimal bet, which is called the "blind". This 
term arises from the fact that they bet before actually 
seeing their cards. The blinds also ensure that there is 
some money in the pot to play for at the very start of 
the game. The blind b increases exponentially with time, 
and typically changes to the value 40 S, 60 $, 100 $, 200 $, 
300$, 400$,... every 10-15 minutes on Internet tourna- 
ments, hence being multiplied by a factor 10 every hour 
or so. We shall see that the growth rate of the blind en- 
tirely controls the pace of a tournament, a phenomenon 
observed in another context in Therefore, the fact 
that the blind grows exponentially with time must be a 
major ingredient of any realistic model of poker. 

• The next players post their bets (> b) according to 
their evaluation of the two cards they each receive. There 
are subsequent rounds of betting following the succes- 
sive draws of five common cards. Ultimately, the betting 
player with the best hand of five cards (selected from 
its two cards and the five common cards) wins the pot. 
Most of the deals end up with a player winning a small 
multiple of the blind. However, during certain deals, two 
or more players can aggressively raise each other, so that 
they finally bet a large fraction, if not all, of their chips. 
This can happen when a player goes all-in, hence betting 
all his chips. Any serious model of poker should take into 
account the fact that players often bet a few blinds, but 
sometimes end up betting all or a large fraction of their 
chips. 

Once a player loses all his chips, he is eliminated. Dur- 
ing the course of the tournament, some players may be 
redistributed to other tables, in order to keep the number 
of tables minimum. 

Retaining the two main ingredients mentioned above, 
we now define a simple version of poker which turns out 
to describe quantitatively the evolution of real poker tour- 
naments. The No initial players are distributed at tables 
with 9 — 10 seats. They receive the same amount of 
chips xq S> bo, where bo is the initial blind. The ratio 
Xo/bo is typically in the range 50 — 100 in actual poker 
tournaments. 

• The players take turns at dealing. In the model, 
only the player next to the dealer, dubbed the "blinder" , 
posts the blind bet. The blind increases exponentially 
with time as, bit) = boexp(t/to). 



• The tables run in parallel. At each table, the players 
receive one card, c, which is a random number uniformly 
distributed between and 1. 

• We define a critical hand value Co <E [0, 1]. The fol- 
lowing players bet the value b with probability e(c), if 
< c < Co. e(c) is an evaluation function, whose de- 
tails will be immaterial. Intuitively, e(c) should be an 
increasing function of c, implying that a player will more 
often play good hands than bad ones. We tried several 
forms of e(c), obtaining the same results. In our simu- 
lations, we choose e(c) = c™, where n is the number of 
players having already bet, including the blinder. In this 
case, e(c) is simply the probability that c is the best card 
among n + 1 random cards. This reflects the fact that a 
player should be careful when playing bad hands if many 
players have already bet. Determining the optimal eval- 
uation function for a given 9, in the spirit of Borel's and 
von Neumann's analysis for 9 — 2, is a formidable task 
which is left for a future study (l(J. 

• The first player with a card c > Co goes all-in, so 
that q — 1 — Co is the probability to go all-in. The next 
players including the blinder can follow if their card is 
greater than Co, and fold otherwise. If a player with a 
card c > Co cannot match the amount of chips of the first 
player all-in, he simply bets all his chips, but can only 
expect to win this amount from each of the other players 
going all-in. 

• Finally, the betting player with the highest card wins 
the pot and the blinder gets the blind back if nobody else 
bets. The players left with no chips are eliminated, and 
after each deal, certain players may be redistributed to 
other tables, in a process ensuring that the number of 
tables remains minimum at all times and that no table 
has less than [N/6] players, where [•] denotes the integer 
part function. After a deal is completed at all tables, 
time is updated to t + 1, and the next deal starts. This 
process is repeated until only one player is left. 

III. POKER MODEL WITHOUT ALL-IN 
PROCESSES 

Let us first consider the unrealistic case q = 0. The 
amount of chips or stack x(t) of a given player evolves 
according to x(t + 1) = x(t)+e(t)b(t). The effective noise 
e(t) should have zero average since all players are consid- 
ered equal and there is therefore no individual winning 
strategy in the mathematical sense. e(t) is also Marko- 
vian, since successive deals are uncorrelated. We define 
x(t) as the statistical average of x(t). If the typical value 
of x ~ x(t) remains significantly bigger than the blind 
b(t), we can adopt a continuous time approach. Hence, 
the evolution of x(t) is that of a generalized Brownian 
walker: 

^=o-b(t)rj(t), (1) 

where a 2 = e 2 is a constant of order unity, and r)(t) is a 6- 
correlated white noise. The number of surviving players 
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with x chips, P(x,t), evolves according to the Fokker- 
Planck equation 

dP _ a 2 b 2 (t) d 2 P 

~m ~ 2 a?' 

with the absorbing boundary condition P(x — 0, t) — 
0, and initial condition P(x,t — 0) = S(x — x ). This 
kind of problem arises naturally in physics in the context 
of persistence, which is the probability that a random 
process x(t) never falls below a certain level [[y} [l2[ [L3| . 
Defining 

a 2 b\to / 2i 

r(t) = — ^ (e*o 



1 



(3) 



Eq. (0) can be solved by the method of images 
No 



P(x,t) = 



2t(*) 



" 2r(i) 



(4) 



V27rr(t) 

Note that the above result holds for any form of b(t), 
provided that one properly defines r(t) = a 2 J Q b 2 {t') dt' . 
For large time (or r ^> 1), the distribution of chips 



(or r > 
becomes scale invariant 



P(x,t) 



Nit) f x 
x(t) 1 \x(t) 
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where the density of surviving players is given by 

N(t) 2x * 
= — = e *o . 

No y/7Tto^(Tb 

We find that the decay rate of the number of players is 
exactly given by the growth rate of the blind, which thus 
controls the pace of the tournament. The total duration 
of a tournament tf is typically 



tf 
to 



1 



ln(A )--ln(t )+ln 



(7) 



which only grows logarithmically with the number of 
players and the ratio xo/bo- The average stack is pro- 
portional to the blind 



x(t) 



Np 
N(t) 



x 



-bit). 



(8) 



When to ^S> 1, this expression implies that x(t)/b(t) » 1, 
hence validating the use of a continuous time approach. 
Finally, we find that the normalized distribution of chips 
is given by the Wigner distribution 

fiX) = ^-Xe^ x \ F(X) = l-e-^ x \ (9) 



where F(X) = j Q f(Y) dY. Equivalently, in the context 
of persistence, / is naturally found to be the first excited 
eigenstate of the quantum harmonic oscillator |l2| . The 
scaling function / is universal, i.e. independent of all the 
microscopic parameters (bo, to, xq...). In Fig. [j], we plot 
the normalized distribution f(X) — x(t)P(x, t)/N(t) and 
its cumulative sum FiX) as a function of X — x/x(t), 
as obtained from extensive numerical simulations of the 
present poker model with q = 0. We find a perfect data 
collapse on the analytical result of Eq. (^|) . 
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FIG. 1: We plot the normalized distribution of chips f(X) 
and its cumulative sum F(X) obtained from numerical simu- 
lations of our poker model with q — (thin lines; iVo = 10000, 
t = 2000, xo/bo = 100, 10000 "tournaments" played). These 
distributions are extracted at times for which N(t)/No = 
50%, 30%, 10%. The dashed lines correspond to the exact 
solution, Eq. (^). 



IV. POKER MODEL INCLUDING ALL-IN 
PROCESSES 

Let us now consider the more realistic case q > (or 
Co < 1). A priori, it seems that q is a new parameter 
whose precise value could dramatically affect the dynam- 
ics of the game. In reality, q must be intimately related 
to the decay rate t^ 1 of the number of players, which 
is imposed by the exponential growth of the blind. To 
see this, let us first compute the decay rate due to the 
all-in processes. At a given table, and for small q, the 
probability that an all-in process occurs is given by 



Pea-t 



q 2 0i6-l) 



(10) 



where the factor q 2 is the probability that two players go 
all-in, and 8(6 — l)/2 is the number of such pairs. Ex- 
pecting q -C 1, we have neglected all-in processes involv- 
ing more than two players. During a two-player all-in 
process, there is a probability 1/2 that the losing player 
is the one with the smallest number of chips (he is then 
eliminated). Cumulating the results of the N/8 tables, 
we find the density decay rate due to all-in processes 

dN 1 N N 

IT = -o X ^- xP all-i„ = -7 , (11) 

Ut all-in I a tail-in 

4 

iall-in = -TTTj 7T- (12) 



q 2 (8-iy 



We now make the claim that the physically optimal 
choice for t a u-m; an d hence for q, is such that the decay 
rate due to all-in processes is equal to the one caused 
by the stack fluctuations of order b(t). Since the to- 
tal decay rate will be shown to remain equal to t$ , 
^all-in = 2to should hold, since inverse decay rates add 
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up. If iaii-in < 2^o, the game is dominated by all-in pro- 
cesses and x(t) can get rapidly large compared to b(t). 
The first player to go all-in is acting foolishly and takes 
the risk of being eliminated just to win the negligible 
blind. Inversely, if £ a ii-in > 2io, players (especially those 
with a declining stack) would be foolish not to make the 
most of the opportunity to double their chips by going 
all-in. We expect that real poker players would, on aver- 
age, self-adjust their q to its optimal value. Finally, we 
find that q is not a free parameter, but should take the 
physical value 



l)to' 



(13) 



We now write the exact evolution equation for the 
number of surviving players with x chips, combining the 
effect of pots of order b and all-in processes 



dP 
~dt 



a 2 b 2 d 2 P 



-(K(P)-P), 



2 dx 2 ' t 
where the non linear all-in kernel K is given by 



(14) 



K(P) 



/2 



N 



P{x-y)^-dy 



1 P(v) 



(15) 



and where we have dropped the time variable argument 
for clarity. In Eq. (|l 
simply the rate of a 



the factor 2/t = q 2 (9 - 1) is 
1-in processes involving the consid- 
ered player, without presuming the outcome of the event. 
In addition, the first term of Eq. ( [l5| ) describes processes 
where the considered player has doubled his chips by win- 
ning against a player with more chips than him. The 
second term corresponds to an all-in process where the 
player has won against a player with less chips than him 
(and has eliminated this player). Finally, the last term 
describes the loss against a player with less chips than 
him (otherwise the considered player is eliminated). In- 
tegrating Eq. (p~5|) over x, we check that the probability to 



survive an all-in process is |, the two first terms adding 
up to | . Indeed, the player survives if he wins (with prob- 
ability ~) or if he loses, but only against a player with 
less chips (with probability j). We recover the decay rate 
associated to all-in processes, (l — |) = t~ n _ in . 

We now look for a scaling solution of Eq 
form 



5h of the 



P(x,t) 



X 



m 



/{ x(t) 



(16) 



where the integral of / is normalized to 1, so that N(t) = 
\/x(t). Plugging this ansatz into Eq. (f5), we find that 




FIG. 2: We plot the normalized distribution of chips f(X) 
and its cumulative sum F(X) obtained from numerical sim- 
ulations of our poker model (thin dotted lines, iVo = 10000, 
t = 2000, xo/bo = 100, 10000 "tournaments" played). These 
distributions are extracted at times for which N(t)/No = 
50%, 30%, 10%. The dashed lines correspond to the numerical 
solution of the exact Eq. ([L8|). The data recorded from 20 real 
poker tournaments (totalizing 1584 players still in) are also 
plotted (full lines), and are found to agree remarkably with 
the present theory. Note that f(X) for real tournaments was 
obtained by differentiating a fitting function to the actual cu- 
mulative sum. We also plot the standard but noisier bin plot 
of the distribution of chips in real poker tournaments (circles). 



one must have x(t) ~ b(t) for all the terms to scale in the 
same manner. Defining 



x{t) 



t ab(t) 



x(t) 



e*o 



(17) 



and the scaling variable X = x/x(t), we obtain the fol- 
lowing integrodifferential equation for f(X) 

f"(X) + Xf'(X) + I/(X/2) /+~ f(Y) dY 
+ J X/2 f(X-Y)f(Y)dY 



+ U + °°f(X + Y)f(Y)dY = 0, 



(18) 



with the boundary condition /(0) = 0. We did not suc- 
ceed in solving this equation analytically. However, the 
small and large X behavior of f(X) can be extracted 
from Eq. 



/(*) - ~, 



X~*+oo 



2[ie 



-fiX 



(19) 



Thus, when including all-in processes, the universal scal- 
ing distribution decays more slowly than for q = 0. 
Eq. ( [l8| ) can be easily solved numerically using a stan- 
dard iteration scheme, and we find /i ~ 1.562. 

In Fig. H we plot the normalized distribution f(X) as 
a function of X = x/x(t) obtained from extensive numer- 
ical simulations of the present poker model, with q given 
by Eq. (p^|). We find a perfect data collapse on the nu- 
merical solution of the exact scaling equation Eq. (p^). 
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FIG. 3: We plot the integrated density distribution F(X) 
extracted from Internet tournaments (same as Fig. [| black 
line), and from the four main events of the WPT 2006 sea- 
son (ranking after day 1 and day 2 totalizing 1256 players; 
turquoise line), and compare them to the numerical solution 
of the exact Eq. ( |l8| ) (dashed line), and to the analytical re- 
sult Eq. (^) of the poker model for q = (dotted line). The 
bottom insert illustrates the exponential decay of 1 — F{x), 
instead of the Gaussian decay predicted by the simple model 
with q = 0. In the top insert, we plot F(X) for each individual 
WPT tournaments (day 1 only; for the sake of clarity). 



In order to check the relevance of this parameter-free 
distribution to real poker tournaments, we visited two 
popular on-line poker playing zones, and followed 20 no- 
limit Texas hold'em tournaments with an initial number 
of players in the range 250 — 800. When the number of 
players was down to the range N ~ 60 — 130, we manu- 
ally recorded their number of chips Fig. || shows the 
remarkable agreement between these data and the results 
of the present model. The maximum of the distribution 
corresponds to players holding around 55% of the average 
number of chips per player. In addition, a player owning 
twice the average stack per player {X = 2) precedes 90% 
of the other players, whereas a player with half the aver- 
age stack (X = 1/2) precedes only 25% of the other play- 
ers. In Fig. [| we compare these results to data collected 
from the four main events of the World Poker Tour 2006 
season Jl8| . Although the level of play is incomparably 
better than on typical Internet poker rooms (the buy-in 
of 10000$ or more is also incomparable), the stacks dis- 
tributions are very similar and decay exponentially (see 
the prediction of Eq. (|l9|)), as illustrated in the bottom 
insert of Fig. H. The model without all-in events (q = 0) 
would predict a faster Gaussian decay. The fact that we 
find similar results for two very different kinds of poker 
tournaments certainly justifies the universal nature of the 
present theory. 
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FIG. 4: We plot the average number of chip leaders Lm 
as a function of the number of initial players No, finding a 
convincing logarithmic growth (full symbols correspond to the 
case q = 0). The insert shows the logarithmic growth of r ma x 
(defined in the text). The dashed lines correspond to log- 
linear fits of the data. 
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FIG. 5: We plot the normalized cumulative distribution G(z) 
of the fortune of the chip leader recorded when N(t)/No = 
40%, 20%, 10% (thin dotted lines, all-in processes included, 
iVo = 10000, t = 2000, xo/bo = 100, 10000 "tournaments" 
played). We also plot the derivative g(z) of a multi- variable 
fit of these data (full line) as well as the standard scatter 
bin plot (circles). The data convincingly follow the universal 
Gumbel distribution (dashed lines). 



V. PROPERTIES OF THE CHIP LEADER 

We now consider the statistical properties of the player 
with the largest amount of chips at a given time, dubbed 
the chip leader. First, we consider the average number 
of chip leaders Ln in a tournament with A^o initial play- 
ers. In many competitive situations |3| 4, 5], arising for 
instance in biological evolution models [3, 4], it is found 
that LAr grows logarithmically with the number of com- 
peting agents TVo, a general result which has been estab- 
lished analytically in H . We confirm that in the present 
model, with or without all-in processes, the same phe- 
nomenon is observed (see Fig. 0). We have also computed 
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the average maximum ratio r max = sup t ahead/a; ■ In the 
present model, Xi eac j/x increases rapidly on a scale of or- 
der to, and then decays (almost linearly with time) to 
~ 1.5, where it becomes non self- averaging due to large 
fluctuations at the end of the tournament. Fig. illus- 



trates the logarithmic growth of 



as a function of 



A^o- For Nq = 500, which is typical of Internet tourna- 
ments, we find r max « 4.6, which is fully compatible with 
a superficial analysis of real data. 

Extreme value statistics have recently attracted a lot 
of attention from physicists in various contexts p4| . 
In this regard, we have checked that z — (xi ea d ~ 
ahcad)/(2; 2 iead — xf cad y^ 2 is distributed according to the 
universal Gumbel distribution 



V6 



exp[— Z — exp(— Z)], 



(20) 



where Z = 7rz/V6 -I- 7, and 7 is Euler's constant. Such 
a behavior, which is typical of independent, or at least 
weakly correlated random variables ji~5| , is illustrated on 

Fig. r 



VI. CONCLUSION 

In this paper, we have developed a quantitative theory 
of poker tournaments and made the connection between 
this problem and persistence in physics, the leader prob- 
lem in evolutionary biology, and extreme value statistics. 
In particular, we have identified the two main ingredients 
controlling the dynamics of a tournament: the exponen- 
tial increase of the blind, and the necessity to include 
all-in events where at least two players bet their entire 
stack. In order to mimic the play of "intelligent" play- 
ers, we found that the probability of going all-in should 
take a well-defined value. This theory leads to a quan- 
titative understanding of the scale-invariant stack dis- 
tribution observed in Internet and WPT tournaments, 
and predicts rich statistical features concerning the chip 
leader. 



In a future work (T(|, we plan to implement in our 
model the optimal strategies for folding, betting or go- 
ing all-in, hence eliminating the only free parameter q. 
Preliminary results fl(| indicate that the optimal prob- 
ability qo to be the first to go all-in is a simple function 
of the current pot P, of the chip stack x of the con- 
sidered player, and of Xk, the stack of the fc-th player 
left to bet (among a total of n such players). Defin- 
ing Xk = min(xfe,x), and q% as the probability that the 
player k calls the all-in bet of the first player (and ne- 
glecting multiple calls), we find O] 



qo 



P = 



x k + P 

X k + 2P 7 



a n (1 - qjY 



(21) 



(22) 




fc=i j=k 

where go is the solution of the implicit Eq. (1221) 
inserting the expression of q^ obtained in Eq. (BIT) 
tailed analysis of Eqs. ( pl| , |22] ) reveals that the obtained 
optimal strategy perfectly reproduces qualitative features 
observed in real tournaments, notably the fact that play- 
ers with a small stack go more often all-in than others 
(and are often called). In addition, direct confrontations 
between two players owning a big stack (in units of x) are 
rare, except if the pot is already huge, and only happens 
when both players have a very good hand. 

Finally, it would be interesting to obtain access to 
the full dynamical evolution of a large sample of real- 
life poker tournaments, in order to check the predictions 
of the model concerning the chip leader and to identify 
other remarkable statistical properties of poker tourna- 
ments. 
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